from numpy import pi, power, exp

def eta(t):
    q = exp(2*pi*(0.0+1.0j)*t)
    res = exp((0.0+1.0j)*pi*t/12.0)
    for n in range(1000):
        res = res * (1-power(q,n+1))
    return res

